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ABSTRACT 

When neutron stars reside in transient X-ray binaries, their crustal layers become heated dur¬ 
ing accretion outbursts and subsequently cool in quiescence. Observing and modelling this 
thermal response has yielded valuable insight into the physics of neutron star crusts. However, 
one unresolved problem is the evidence in several neutron stars for an extra energy source, 
located at shallow depth in the crust, that is not accounted for by standard heating models. 
Its origin remains puzzling, and it is currently unclear whether this additional heating occurs 
in all neutron stars, and if the magnitude is always the same. Here, we report on Chandra 
observations that cover two years after the 2012 outburst of the transient neutron star X-ray 
binary Swift J174805.3-244637 in the globular cluster Terzan 5. The temperature of the neu¬ 
tron star was elevated during the first two months following its ~8 week accretion episode, 
but had decayed to the pre-outburst level within ~100 d. Interpreting this as rapid cooling of 
the accretion-heated crust, we model the observed temperature curve with a thermal evolution 
code. We find that there is no need to invoke shallow heating for this neutron star, although 
an extra energy release up to ~1.4 MeV nucleon”^ is allowed by the current data (2a confi¬ 
dence). We also present two new data points on the crust-cooling curve of the 11-Hz X-ray 
pulsar IGR J17480-2446 in Terzan 5, which was active in 2010. The temperature of this neu¬ 
tron star remains significantly above the pre-outburst level, but we detect no change in the 
thermal emission since the previous measurements of 2013 February. This is consistent with 
the slow crust cooling expected several years post-outburst. 

Key words: stars: individual (Swift J174805.3-244637) - stars: pulsars: individual (IGR 
J17480-2446) - stars: neutron - X-rays: binaries - globular clusters: individual: Terzan 5 


1 INTRODUCTION 
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Low-mass X-ray binaries (LMXBs) are excellent laboratories 
to gain ins ight into the structure an d composition of neutron 
stars (e.g., iMilleil I 2 OI 3 I: lOzell 1201 1 for recent reviews). One 
promising avenue is to study how neutron stars are heated dur¬ 
ing transient accretion events and cool thereafter. This ther- 
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Table 1. Overview of neutron star LMXBs for which crustal cooling has been monitored. 


Source name 

*0 

(MJD) 

^ob 

(yr) 

i^7ob 

(gs-l) 

D [ref] 
(kpc) 

^spin 

(Hz) 

Povh [ref] 

(hr) 

kT^ 

base 

(eV) 

Remarks 

Cooling 

references 

KS 1731-260 

51930.5 

12.3 

1.9 X 10^'' 

7.0 [1] 

524 [2] 



Superburst source 

3,4,5,6 

MXB 1659-29 

52159.5 

2.4 

1.1 X lOi'^ 

10 [7] 

567 [8] 

7.1 [9] 


Eclipsing 

5,10,11,12,13 

XTE J1701M62 

54322.0 

1.6 

1.1 X lOi® 

8.8 [14] 




Transient Z-source 

15,16 

EXO 0748-676 

54714.0 

24 

4.8 X 10^® 

7.1 [17] 

552 [18] 

3.8 [19] 

94.6 

Eclipsing 

20,21,22,23 

IGR J17480-2446 

55556.0 

0.17 

2.3 X lO^'^ 

5.5 [24] 

11 [25] 

21.3 [25] 

73.6 

In Terzan 5 

Transient Z-source 
X-ray pulsar 

26,27,28,29 

MAXI J0556-332 

56052.1 

1.3 

9.1 X lOi’’ 

45 [30] 


16.4/9.8 [31] 


Transient Z-source 

30 

Swift J174805.3-244637 

56166.0 

0.15 

9.6 X lO^® 

5.5 [24] 



89.9 

In Terzan 5 

29 


Note. - The assumed end of the outburst is indicated by Iq, whereas fob ^nd Mob estimates of the outburst duration and average mass-accretion rate, 
respectively. Any available information on the distance (D), neutron star spin frequency (t'spin). orbital period (Porb). and temperature measured before 


MAXI J0556-332. References: 1 =lMuno et alJ i200Gll. 2=ISmith et alJ il997ll. 3 =IWiinanc 

setalJ J2001bli. 4 = 

Wiinands et alJ i2003h. 5 = 

Cackett et alJ 

l200d). 6 =ICackett et alJ J2010bl). 7 =lMuno et alj J200lh, 8 =lwiinands et alj J2001ah. 9 = 

Cominskv & Wood ' 

1984h. 10 = Wiinands et al, 

II 

o 

o 

<N 

Iwiinands et alJ i2004h. 12 =ICackett et al 

i2008li. 13 =ICackett et alJ i2013li. 14 =lLin et al 

J2009h. 15 =lFridriksson et alJ J201(t). 16 =lFridriksson et alJ 

boi lb. 17 =lGallowav et al.l b008t). 18 = 

Gallowav et alj J2010h. 19 =|Parmar et alj Jl986ll. 20 = Deaenaar et alj J2009h. 21 =lDeaenaar et al.| |2011aJ). 22 


IPegenaa^t alJbOllblT 28 = lDegenaar et alj b013bl) . 29 = this work. 30 = lHoman et aljbOl^ . 31 = ICome isse et alj j2012h . 


mal response can give detailed information about the struc¬ 
ture and composition of the crust, as w ell as the density and 


superfluid properties of the core (e.g.. 

Rutledge et alJ l2002bl: 

Shtemin et al. 

2OO7I: Brown & Cummind 

1 20091: IPase & Reddv 

2OI2I: ITurlione 

et al.l 1201.51: IWiinands et al 

201.31: iHorowitz et al. 

2015l:lMedin & Cumming||2015l). 


Transient LMXBs undergo occasional outbursts during which 
matter from the <1 Mq companion is rapidly accreted on to 
the neutron star. The accumulation of matter compresses the 
stellar crust and causes a series of nuclear reactions that in- 
clude electron captures a nd density-driven fu sion reactions (e.g., 
iHaensel & Zdunilall99(l l2008l : ISteineill20l3) . The resulting en¬ 
ergy releas e of MeV per accreted nucleon causes the crust 
to heat up l lBrown et al.lll99d) . creating a temperature profile that 
depends on the location and magnitude of the heat sources, as 
well as the duration and brightness of the accr etion episode (e.g., 
iBrown & Cummin3l2009l : IPage & Reddvll2oT^ . Typical outbursts 
reach a luminosity of Lx — erg (0.5-10 keV) and last 

for weeks to months (e.g., such as seen in Aql X-1, SAX J1808.4— 
3658, and EXO 1745-248), although a handful of neutron stars 
continue to accrete for years to decades at this level (such as e.g., 
MXB 1659-29, EXO 0748-676, and HETE J1900.1-2455). 


Outbursts are usually separated by long (years to decades) 
periods of quiescence with a much lower X-ray luminosity of 
Lx — 10®^“®"^ ergs“^, suggesting that the mass-accretion rate 
on to the neutron star is strongly reduced. The heated crust can 
then cool by thermally conducting the energy gained during the 
outburst into the neutron star core and towards the surface, where it 
is radiated as neutrinos and (X-ray) photons, respectively. This ther¬ 
mal relaxation depends on the structure, thickness, and composition 
of the crust since th e se set the thermal tran s port properties (e.g., 
Lattinjer^ 1 19941: [Rutledge et al.l l2 002bl: iBrown & Cummind 


20091 : Medin & Gumming 20151 : Horowitz et 2015h . as well as 
the temperature of the stellar core, which relate s to its density 
and s uperfluid properties (e.g.. lShternin et alJl2007l : |Page & Reddvl 
l2013t) . 


The X-ray emission observed during outbursts is dominated 
by the accretion flow and hence heating of neutron stars cannot 
be directly observed. However, thermal X-ray emission from the 
incandescent neutron star can often be detected in quiescence. This 


offers the opportunity to probe its thermal response to the accretion 
episode and to monitor the subsequent cooling. 


In the past decade, crust cooling has successfully been moni¬ 
tored for six transiently accreting neutron stars (Table [TJ. In these 
LMXBs, the temperature of the neutron star has been seen to de¬ 
crease systematically for several years following an accretion out¬ 
burst. Comparing these observations with crust-cooling simulations 
has provided valuable new insight into the properties of neutron star 
crusts. For example, it has been established that heat is conducted 
rapidly and hence that t he atomic nuclei in the crust must have a 
highl y ordered structure iShtemin et al.ll2007 ; Brown & Cummin3 


1 20091 : IPage & Reddvl l2()13l : iTurlione et alj 20151) . FurthermoreT 
crust-cooling studies have the potential to place st rong constraints 
on th e composition of the outer c rustal layers (e.g.. lwiinands et aP 
I2OI3I : iMedin & Cummin3 l2015ll . to test the occurrence of non- 
sphe rical nuclear shapes ( ‘nuclear pasta’) in the deep inner crust 
(e.g.. lHorowitz et ^l2015h . and even to probe the superfluid prop- 
erties of the ultradense stellar core (e.g., IPage & Redd^l2013h . 


Despite these successes, several open questions remain. One 
of the current puzzles is that modelling of the crust-cooling curve 
points, in some sources, to the presence of an extra source of en¬ 
ergy (possibly as high as a few MeV nucleon'^) that is not ac¬ 
counted for by current nuclear heating models, but must be located 
at a shallow d epth in the crust (e.g., in MXB 1659-29 and IGR 
J17480-2446: [Brown & CummiM [20091 : IPegenaar et akll^Ol Ibh . 
Understanding the origin of this shallow heating is of particular 
interest because it may also explain the occurrence of highly en¬ 
ergetic thermonuclear bursts called ‘superbursts’, which have igni¬ 
tion properties that require a considerably higher crust temperature 
than can be achieved with curren t nuclear heating models (e.g., in 
4U 1608-52 and EXO 1745-2 48: ICumming et ali2nO^:lKeek et al] 
l2008l : lAitamirano et al.ll201^ . Furthermore, the transition between 
stable and unstable thermonuclear burning on the surface of neu¬ 
tron stars, which gives rises to mHz quasi-periodic oscillations, is 
observed to occur at an accretion rate that is about an order of mag¬ 
nitude lower than predicted by theory (e.g. , in Aql X-1, 4U 1608- 


52, 4U 1636-53, and IGR 117480-2446: Revnivtsev et al. 

I2OOII: 

lAltamirano et alJl2008l:lKeek et alJl2009l:lLinares et al.ll2012 

). This 


discrepancy could also be reconciled if additional crustal heating 
occurs. 
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Terzan 5 


Swift J 174805.3-244637 

^(T5X-3), 


EXO 1745-248 i 
(T5X-1) , 


0 


IGRJ17480-2446 
(T5 X-2) 


2 " 


Figure 1. Accumulated three-colour Chandra!A.C\^ image of the inner 
~1.1 ai'cmin xl.l arcmin of the globular cluster Terzan 5 (0.3-8 keV). 
Red colour coiTesponds to 0.3-1.5 keV, green to 1.5-2.5 keV, and blue to 
2.5-8 keV. This image was constructed using 16 observations in which all 
bright X-ray transients were quiescent, amounting to 620 ks of total expo¬ 
sure time. The positions of the three known transient neutron star LMXBs 
are indicated by circular regions with 1.5 arcsec radii. 


Possible explanations for the origin of extra heating in- 
clude additional el ectron captures at shallow crustal layers (e.g., 
lEstrade et alJ2Ql ih . and previously unacco unted nuclear fusion re - 
actions occurring deep in the crust (e.g., iHorowitz et alj [2008h . 
However, electron captures may be too weak and nuclear fusion 
may occur at too high density to be r econciled with observat ional 
constraints on the extra heat release toeeenaar et alj[2013ah . An 
alternative mechanism that may be consistent with observations is 
a convective heat flux driven by the separation of light and heavy 
nuclei in the outer crustal la yers dMedin & Cumminal201 lll2015l : 
IPegenaar et al.ll^l3aLl2014ll . Currently it is unclear whether extra 
crustal energy release occurs for all neutron stars, and if the depth 
and magnitude are always the same. Performing crust-cooling stud¬ 
ies for additional sources, in particular probing a range of source 
parameters (e.g., outburst duration and brightness, orbital and spin 
period) can potentially shed light on this unresolved problem. 

Crust-cooling studies initially focused on neutron stars ex¬ 
hibiting long (>1 yr) accretion outbursts to ensure that signifi¬ 
cant heating occurred and hence optimizing the chances of de¬ 
tecting the subsequent cooling. However, it has been demonstrated 
that significant heat ing can also occur during shorter outbursts (in 
IGR J17480-2446; IPegenaar & Wiinandsll20lT3 : IPegenaar et al.l 
I 2 OI lbl . [2013bh . Following neutron stars after weeks to months long 
accretion outbursts may be particularly suited to investigate the 
presence of shallow heat, since the su bsequent cooling should be 
faster (e.g., iBrown & Cummin3l2009ll and detailed modelling of 
the complete cooling curve can therefore start sooner. Moreover, 
short outbursts are much more common than prolonged ones and 
hence may offer better opportunities to increase the number of ob¬ 
served crust-cooling curves. 


1.1 Swift J174805.3-244637 in Terzan 5 

Terzan 5 is a strongly absorbe d (A^h — 10^^ cm~^: lHeinke et ^ 
l2006bl : [Bahramian et alj|2014h Milky Way globular cluster that is 
located at an estimated dis tance of ~4.6-8.7 kpc (e.g.. lCohn et al.l 
l2002l : [0rtolani et al]|2007h . It harbours a large number of faint X- 
ray point sour ces, many of which are candidate quiescent neutron 
star LMXBs l lHeinke et al .1120065 . see also Fig. [TJ. Indeed, two 
confirmed transient neutron star LMXBs were previously identi¬ 
fied in the cluster: EXO 1745-248 and the 11-Hz X-ray pulsar IGR 
J17480-2446. The former was active at least in 2000, 2011, and 
2015 (e.g., Heinke et alj20'03; Altamirano et ^J2012|; Barret|20ld: 


ISerino et al. I I 2012 I; Iai 


tamirano et al.l l2015l: iTremou et alj l2015h 

wherea s the latter was responsible for the 2010 outburst of Terza n 
5 (e.g.. IPapitto et ^1201 ll : iMiller et'^l201 ll : iMotta et'n]|201lh . 
X-ray activity was also ob served from the cluster in 19 80, 1984, 
1990, 1991, and 2002 (see IPegenaar & Wiinandsll2012l for a his¬ 
toric overview up until 2012), but due to the high source density 
and lack of sub-arcsecond spatial resolution observations, it cannot 
be pinpointed with certainty which object(s) caused these outbursts 
(see Fig.[TJ. In 2012, a third transient neutron star LMXB was dis¬ 
covered in Terzan 5: Swift 1174805.3-244637 (also known as T5 
X-3: lBahramian ^ alj2014ll . The cluster core with the locations of 
the three transient LMXBs is shown in Fig.[T] 

The 2012 discovery outburst of Swift 1174805.3-244637 
started in early luly and the source remained active for ~7- 
8 weeks. Puring this time it was detected at a mean 0.5-10 keV lu¬ 
minosity of Lx — 9x10^® ( D/5.5 kpc)^ erg s~^, an d it displayed 
a thermonuclear X-ray burst l lBahramian et alj[2014ll . Analysis of 
archival Chandra data revealed that the pre-outburst quiescent X- 
ray spectrum was dominated by thermal emission from the neu¬ 
tron star surface, but with an ~20-30 per cent contribution from a 
non-thermal emission component. The neutron star temperature re¬ 
mained constant between 2003 and 2012, but smal l variations in the 
quies cent non-thermal emission were observed jBahramian et al.l 
|2014|) FI Such a hard (variable) emission tail is often detected for 
quiescent neutron star LMXBs and is typically modelled by a sim¬ 
ple power law with an index of F ~ 1 — 2. It is thought to arise 
from a different emission process, possibly related to the presence 


star (e.g.. Camnana et al. 19981; lonker et al. 2004 

; Wiinands et al. 

2 OO 5 I: ICackett et al.ll2010al; iDegenaar et al.ll2012l; 

Bemardini et al. 

2013l:IChakrabartv et alj|2014l). 


When the 2012 outburst of Swift 1174805.3-244637 ceased 
and the source transitioned to quiescence, we took the opportunity 
to search for crustal cooling in this neutron star. 


2 DATA ANALYSIS AND RESULTS 

2.1 Observations and analysis procedures 

Terzan 5 has been observed five times since the end of the 2012 
outburst of Swift 1174805.3-244637 as part of a Chandra Target 
of Opportunity programme to monitor the thermal evolution of the 
neutron star. Two additional Chandra observations were performed 


^ No outburst was detected from Swift 1174805.3-244637 in the epoch 
2003-2012, s uggesting that all these observations sampled the same quies¬ 
cent episode. iBahramian et al.l i2014h did propose that because of the rel¬ 
atively high temperature of the neutron star. Swift 1174805.3-244637 may 
have been responsible for one or more of the historic outbursts of Terzan 5 
for which no accurate position information was available. 
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Table 2. ChandralACIS-S observations of Terzan 5. 


ObsID 

Date 

Exposure Time 
(ks) 

Count rate 
(10-3 cs-l) 


Pre-outburst observations 


3798 

2003 Jul 13/14 

39.5 

10.0 ± 1.0 

10059 

2009 Jul 14/15 

36.4 

7.3 ± 0.5 

13225 

2011 Feb 17 

29.7 

5.0 ±0.4 

13252 

2011 Apr 29/30 

39.5 

6.9 ± 0.4 

13705* 

2011 Sep 5 

13.9 

5.5 ± 0.6 

14339* 

2011 Sep 8 

34.1 

6.0 ±0.5 

12454 

2011 Nov 3 

9.8 

7.3 ± 0.8 

13706 

2012 May 13 

46.5 

7.7 ±0.4 


Post-outburst observations 


14475 

2012 Sep 17/18 

30.5 

10.0 ±0.6 

14476 

2012 Oct 28 

28.6 

9.4 ±0.6 

14477* 

2013 Feb 5 

28.6 

7.4 ±0.5 

14625* 

2013 Feb 22 

49.2 

7.1 ± 0.4 

15615* 

2013 Feb 24 

84.2 

6.4 ±0.3 

14478 

2013 Jul 16/17 

28.6 

7.3 ± 0.5 

14479* 

2014 Jul 15 

28.6 

5.7 ±0.5 

16638* 

2014 Jul 17 

71.6 

5.6 ± 0.3 

15750* 

2014 Jul 20 

23.0 

6.1 ± 0.5 


Note. - The count rates for Swift J174805.3-244637 are given for the 
0.3—10 keV energy range. Quoted uncertainties are at the Itr level of 
confidence. Spectra extracted from observations spaced close in time 
are marked by an asterisk and were combined to improve the statistics. 


to stud y the 11-Hz X-ray pulsar in quiescence foegenaar et al.l 
l2013bh . Together these seven observations cover ~2 yr after the 
discovery outburst of Swift J174805.3-244637 (Tabled. For all 
observations the ACIS-S3 chip was used with the ‘faint’ timed data 
mode, either in full frame or in a 1/4 sub-array. None of the obser¬ 
vations suffered from background flares. We reduced and analysed 
the data using the CIAO tools version 4.4 and CALDB version 4.5.2. 

Source events were obtained by using a circular region with 
a radius of 1.5 arcsec (see Fig. [T]l. A source-free circular region 
with a radius of 40 arcsec placed ~2 arcmin from the cluster core 
was used to extract background events. For each observation we 
extracted count rates using DMEXTRACT (Table |3, whereas the 
spectra and response files were generated with SPECEXTRACT. Af¬ 
ter grouping the spectra with GRPPHA into bins of >15 photons, 
we fi tted the data i n the 0.3-10 keV range using XSPEC (version 
12 .8 ;|a maud|[l99^ . Through out this work we ass ume a distance 
of D = 5.5 kpc for Terzan 5 dOrtolani et alj|200% . and quote un¬ 
certainties at the la level of confidence. 

During the three exposures that were taken within ~ 20 d from 
one another in 2013 February, the count rates of Swift J174805.3- 
244637 were similar (Table 1^. To improve the statistics we com¬ 
bined these three spectra (and the weighted response files) using 
COMBINE_SPECTRA. Similarly, the three exposures obtained in 
2014 July had comparable count rates and were therefore combined 
(Table |2ll. This resulted in post-outburst spectra for five different 
epochs. 

To study any possible thermal response of the neutron star to 
the 2012 outburst, we also analyse the sp ectra of all pre-outburs t 
observations in this work (Table[2 see also lSahramian et al.l2014h . 
These data were reduced in the same way as detailed above. The 
two exposures of 2011 September (Tablej^ were combined to cre¬ 
ate a single spectrum using COMBINE_SPECTRA. This left us with 
seven pre-outburst spectra. 


2.2 X-ray spectral analysis 

Following [Bahramian et alj ( l2014h . we use d a combination of a 
neutron star atmosphere model iNSATMOS: iHeinke et al.ll2006ah . 
and a power-law (PEGPWRLW, where we set the normalization to 
represent the unabsorbed 0.5-10 keV flux) to fit all spectra. For the 
NSATMOS component we fix the neutron star mass (M = 1.4 Mq) 
and radius (R — 10 km), the distance (D = 5.5 kpc), and the nor¬ 
malization (Ansatmos = 1, i-c., the emission radius was assumed 
to correspond to the entire neutron star surface). The data did not 
provide sufficient statistics to constrain variations of the power-law 
index when it was allowed to vary between the different observa¬ 
tions. We therefore assumed this parameter to be constant over the 
entire data set and only allowed the power-law normalization to 
vary. To account for interstellar absorption we included the TBABS 
model, adopting the VERN cross-sec tions and WILM abundances 
dVerner et aljl99^ : lwilms et al.l200^ . The hydrogen column den¬ 
sity was assumed to remain constant at all epochs (i.e., this param¬ 
eter was tied between all data sets), which seems reasonable since 
the source did not display any eclipses o r dips and is therefo re likely 
not viewed at very high inclination (see lMiller et al.ll2009l) . 

We first addressed whether there was evidence for variability 
among the five post-outburst spectra. Assuming that both the ther¬ 
mal and the non-thermal model components were the same at all 
epochs resulted in a poor fit (xi = 1-51 for 105 dof; with a p-value 
of = 5.9 X 10““*). Allowing the normalization of the power 
law to vary improved the fit (x2 = 1-20 for 101 dof, = 0.08). 
However, a more significant improvement was achieved by allow¬ 
ing the temperature to vary (with the power-law normalization tied; 
xi = 1.10 for 101 dof, Px ~ 0.22). Leaving the power-law nor¬ 
malization free in addition to the temperature, further improved the 
fit iXv ~ 1-02 for 97 dof, = 0.41), although an P-test suggests 
a probability of 0.02 that such an improvement occurs by chance. 
In Fig. 12 we show the X-ray spectra of the first (2012 September) 
and last (2014 July) post-outburst observations. This plot illustrates 
that the largest spectral changes occurred below ~3 keV, where the 
thermal emission from the neutron star dominates (dashed curves 
in Fig. H. We conclude that there is strong evidence for a changing 
neutron star temperature after the 2012 outburst, but there were no 
large variations in the non-thermal power-law emission. 

We proceeded to fit the entire set of twelve (seven pre¬ 
outburst, five post-outburst) quiescent observation s, each with their 
own r esponse files, together in XSPEC. Since [Bahramian et al.l 
( l2014h found no evidence for changing thermal emission before the 
outburst, we only allowed the neutron star temperature to vary inde¬ 
pendently for the new data (i.e., it was assumed to be the same for 
all pre-outburst observations). We determined for each observation; 
the neutron star temperature (as observed at infinity A:T°°)0the to¬ 
tal unabsorbed model flux Px, the power-law flux Px.pow, and the 
relative contribution of the power-law flux /pow (all 0.5-10 keV), 
as well as the thermal bolometric flux Pth.boi (0.01-100 keV). The 
fluxes and corresponding errors were calculated using the CFLUX 
command in XSPEC. The results of our spectral analysis are sum¬ 
marized in Table|2 and the flux evolution is shown in Fig.|2 

Our approach resulted in a good fit (xt = 1.03 for 162 
dof, P^ = 0.37). We obtained a hydrogen column density of 

^ The temperature seen by an observer at infinity relates to the fitted tem¬ 
perature, kT, as A:r°° = kT/{l + z), where 1 + 2 = (1 — Rs/R)~^^^ = 
1.31 is the gravitational redshift for our choice of M = 1.4 Mq and 
R = 10 km. Here Ra = 2GM/c^ is the Schwarzschild radius, G the 
gravitational constant, and c the speed of light. 
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Table 3. Spectral analysis results for Swift J174805.3-244637 in quiescence. 


Epoch 

MJD 

kT°° 

(eV) 

O 

rX 

-^X,pow 

erg cm~^ s“ 

^ th,bol 

/pow 

Lx 

(10^3 

-^th,bol 

erg s~^) 

2003 Jul 

52833.5 

89.7 ± 1.7 

5.0 ±0.6 

Pre-outburst 

1.53 ±0.17 

3.7 ±0.8 

0.30 ±0.16 

1.8 ± 

0.2 

1.4 ±0.3 

2009 Jul 

55027.5 

89.7 ± 1.7 

3.7 ±0.5 

1.02 ± 0.14 

3.7 ±0.8 

0.28 ±0.19 

1.3 ± 

0.2 

1.4 ±0.3 

2011 Feb 

55609 

89.7 ± 1.7 

3.2 ±0.4 

0.58 ± 0.12 

3.7 ±0.8 

0.18 ±0.25 

1.2 ± 

0.2 

1.4 ±0.3 

2011 Apr 

55680.5 

89.7 ± 1.7 

3.7 ±0.4 

0.89 ± 0.13 

3.7 ±0.8 

0.24 ±0.19 

1.4 ± 

O.l 

1.4 ±0.3 

2011 Sep 

55810.5 

89.7 ± 1.7 

2.9 ±0.4 

0.61 ± O.ll 

3.7 ±0.8 

0.21 ±0.22 

1.0± 

O.l 

1.4 ±0.3 

2011 Nov 

55868 

89.7 ± 1.7 

3.8± l.O 

0.78 ±0.30 

3.7 ±0.8 

0.21 ±0.48 

1.4± 

0.4 

1.4 ±0.3 

2012 May 

56060 

89.7 ± 1.7 

4.1 ±0.5 

1.22 ± 0.14 

3.7 ±0.8 

0.30 ±0.16 

1.5 ± 

0.2 

1.4 ±0.3 

2012 Sep 

56187.5 

98.2 ± 2.2 

4.9 ±0.6 

Post-outburst 

0.86 ±0.08 

5.1 ±0.7 

0.18 ±0.23 

1.8 ± 

0.2 

1.8 ±0.2 

2012 Oct 

56228 

93.9 ± 2.6 

4.5 ±0.6 

1.18±0.20 

4.3 ±0.7 

0.26 ±0.21 

1.6 ± 

0.2 

1.5 ±0.3 

2013 Feb 

56340 

90.0 ± 1.6 

3.5 ±0.3 

0.72 ± 0.08 

3.6 ±0.4 

0.21 ±0.14 

1.3 ± 

O.l 

1.3 ±0.2 

2013 Jul 

56489.5 

90.1 ± 2.7 

4.1 ±0.6 

1.18 ± 0.19 

3.8 ±0.7 

0.29 ±0.21 

1.5 ± 

0.2 

1.4 ±0.2 

2014 Jul 

56855.5 

88.0 ± 1.9 

3.3 ±0.3 

0.78 ±0.09 

3.3 ±0.4 

0.24 ±0.15 

l.l± 

O.l 

1.2 ±0.2 


Note. - Tx and Tx.pow represent the total unabsorbed model flux and the power-law flux in the 0.5-10 keV band, respectively. The parameter 
/pow represents the fractional contribution of the power-law component to the total unabsorbed 0.5-10 keV flux. Ftu i,ol the flux of the thermal 
component in the 0.01-100 keV range. Lx and Ltn bol are the total 0.5-10 keV and thermal 0.01-100 keV luminosities for D = 5.5 kpc, 
respectively. The following parameters were kept fixed during the fits: M = 1.4 Mq, = 10 km, D = 5.5 kpc, and A^nsatmos = 1. The 
simultaneous fit resulted in Nn = (2.63 ± 0.10) X 10^^ cm“^, T = 1.7 ± 0.2, and xi = 1-03 (for 162 dof). 



12 5 10 

Energy (keV) 


Figure 2. Comparison between the post-outburst X-ray spectra of 2012 
September (black, filled circles) and 2014 July (red, open squares) for 
Swift J174805.3-244637. The solid lines represent fits using a model con¬ 
sisting of a neutron star atmosphere (dashed curves) and power law (dotted 
curves). The bottom plot shows the data to model ratio of both spectra. 




Time since end of accretion outburst (days) 


A^h = (2.63 ± 0.10) X 10^^ cm“^, and a power-law index of 
r = 1.7di0.2. Both are consistent with the pre-outburst analysis of 
iBahramian et alj j2014l) . All quiescent spectra of Swift J 174805.3- 
244637 are dominated by thermal emission, which contributes 
~ 70 — 80 per cent to the total unabsorbed 0.5-10 keV flux. The 
fractional contribution of the power-law component does not show 
large variations between the different observations (see Fig. 0 and 
Table|^. In fact, this parameter can be adequately fitted with a con¬ 
stant function, yielding /pow = 0.24 ± 0.01 in 2012-2014. 

For the thermal emission component we measure a pre¬ 
outburst temperature of = 89.7 ± 1.7 eV. During the 

first two new observations, obtained ~2 weeks and ~2 months 
after the 2012 outburst, the temperature was elevated to kT°° = 
98.2 ± 2.2 eV and kT°° = 93.9 ± 2.6 eV, respectively. The subse¬ 
quent three observations, obtained ~0.5-2 yr after the outburst, are 
all consistent with the pre-outburst temperature (Table [3}. Indeed, 


Figure 3. Flux evolution in the 0.5-10 keV band after the 2012 outburst 
of Swift J174805.3-244637. From top to bottom: the total unabsorbed flux 
(Lx), thermal flux (Lx,th). power-law flux (Lx, pow). and power-law frac¬ 
tion (/pow). The vertical dotted lines and the numbering on top indicate the 
different years of the monitoring campaign. 

fixing the temperatures of the last three observations to those of the 
pre-outburst spectra provides an acceptable fit (x^ = 1.03 for 165 
dof, = 0.37). 

An L-test suggests that leaving the temperature of the second 
observation free yields a ~ 1.5 ct improvement compared to fixing 
it to the pre-outburst spectra. Leaving the temperature of the first 
observation free results in an ~ 4a improvement. We conclude 
that there is firm evidence for an enhanced temperature during the 
first post-outburst observation, but only a marginal indication that 
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Figure 4. Thermal evolution of Swift J174805.3-244637 compared to 
model calculations. The observations can be described with standard 
physics input (black solid curve), but a shallow heat source up to 
1.4 MeV nucleon”^ is allowed at an approximate 2cr level (dashed curve). 
The dotted curve shows a model calculation with a lower base level, in 
which case a shallow heat source of 0.85 MeV nucleon”^ is preferred by 
the data (see Section lT^ for details). The grey horizontal bar represents the 
temperature measured prior to the 2012 outburst. The vertical dotted lines 
and the numbering on top indicate the different years in the campaign. 


it was still enhanced during the second observation. Nevertheless, 
Fig.|4]illustrates that there is a systematic decrease in temperature 
following the 2012 outburst. Such a temperature evo lution is not 
seen in the pre-outburst data dBahramian et al.lf2Q14h . In analogy 
with that seen in other sources (see Section [231 >. this may point to 
crust cooling. It would imply that the neutron star crust was heated 
during the ~ 8 week accretion outburst but cooled rapidly (within 
~100 d) afterwards. 

In Section we compare the observed temperature curve 
of Swift J174805.3-244637 with thermal evolution simulations to 
study the crust properties of this neutron star. In Section 12.41 we 
also report on empirical decay fits to the temperature curve to pro¬ 
vide a comparison with other crust-cooling sources. 


2.3 Crust-cooling simulations 

To test the hypothesis that the crust of Swift J174805.3-244637 
was significantly heated during its 2012 outburst and that our 
Chandra observations track the subsequent cooling, we performed 
simulations with the neu tron star thermal evolution code of 
[Brown & Cummin3 (l2009l) . using all the physics input and mod¬ 
elling approach described in that work[3 In brie f: the crust compo¬ 
sition was assumed to match the calculations of lHaensel & Zduiii^ 
lll990l) . and the temperature in the atmosphere during outburst was 


® This code is available at: https://github.com/nworbde/dStar. 


set to Tob = 2.8 X 10® K. This is appropriate if H/He burning oc¬ 
curs, as indicated by the propert ies of a type-I X-ray bu rst observed 
from Swift J174805.3-244637 ( iBahramian et al.ll2014h . 

The first and the last detection of the 2012 outburst (with 
Swi///XRT) were on July 6 and August 24, respectively. The source 
was not seen active yet on Ju ne 30, and on August 30 it had 
faded into the local background jBahramian et al.f2014h . The 2012 
outburst thus had a duration of ~49-61 d (fob =0.15±0.02 yr), 
and ended between 2012 August 24 and 30. During this time 
it was detected at an average 0.5-10 keV luminosity of Lx — 
9 X 10®® (D/5. 5 kpc)^ erg s~® . Assu ming a bolometric correc¬ 
tion factor of 2 din T Zand et al.l 120070 . we estimate an average 
mass-accretion rate during the outburst of Mob = RL/GM ~ 
1.5 X 10“® Mq yr“® ~ 9.6 x 10®® g s“® (for M = 1.4 M© and 
D = 10 km). Fitting the 0.5-10 keV luminos ity light curve of the 
outburst presented in IBahramian et al.l (12014l) to a simple broken 
linear decay function suggests that the source faded into the XRT 
background ~36.5 d after the outburst peak, which was observed 
on MID 56129.1. We therefore tentatively set the outburst end to 
to = MID 56166. 

The temperature of the neutron star core was set to Toots = 
1.55 X 10® K to match the observed pre-outburst surface tem¬ 
perature (kT°° =89.7 eV). We assumed an impurity param¬ 
eter Qimp = 1, which corresponds to a highly conductive 

crust as appropriate for all other crust-cooling neutron stars (e.g., 
[Brown & Cumming 2009: ^ge & Redd^ 1201 3l : IPegenaar et al.l 
201 Ihl . |2014 iHorowitz et aiJ l2015h . However, we note that as 


a result of the short outburst duration the cooling curves of 
Swift J174805.3-244637 are not sensitive to this parameter. This 
is because only the outer layers are heated, where the thermal con¬ 
ductivity is govern ed by electron-ion rathe r than electron-impurity 
scattering (see also lPegenaar et aDl201 Ibl) . 

Fig. a shows that the observed temperature curve of 
Swift J174805.3-244637 can be reproduced using the standard 
physics input (solid curve), and does not require the inclusion of 
an additional source of shallow heat (Tablea. We briefly explored 
to what extent the presence of such shallow heating is allowed by 
the current data, using a grid search technique on (Jshaiiow For our 
standard model (i.e., no shallow heating), we found = 0.4 for 
4 dof. Adding any source of shallow heat increases the value, 
hence gives a worse fit. By increasing the magnitude of the shal¬ 
low heat until we obtain a fit statistic of Ax^ — 4 compared to 
the standard model, we can set an approximate 2a upper limit of 
Qshaiiow < 1-4 MeV nucleon“® (xl. = 1.5 for 4 dof; Table |4]l. 
Fig. a illustrates that shallow heating most prominently increases 
the temperature at early times after the outburst (dashed curve). It 
is therefore that observations within ~ 2 months after the end of 
an outbu rst are instrumental to constrain shallow crust he ating (see 
also e.g.. lBrown & Cumming||200^ : [Page & ReddvlEoi2l) . 

We probed the effect of uncertainties in the outburst prop¬ 
erties of Swift J174805.3-244637 on our conclusions about the 
shallow heating. Allowing for any realistic longer/shorter outburst 
duration (fob = 0.13 — 0.17 yr), a factor 2 higher/lower mass- 
accretion rate, or shifting the outburst end 3 d forward or back¬ 
ward (to = MID 56163-56169), changes the obtained limit on 
Qshaiiow by ±0.3 MeV nucleon~® at most. However, a stronger 
effect is caused by changing the assumed base level. For instance, 
lowering the core temperature to Toots = 1.45 X 10® (corre¬ 
sponding to an observed surface temperature of ~87 eV) actu¬ 
ally prefers the addition of shallow heating, yielding a best fit 
for Qshaiiow = 0.85 MeV nucleon”® (x^ = 0.1 for 4 dof; Ta¬ 
ble |4]l. This model is shown as the red dotted curve in Fig. |4] 
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Table 4. Model parameters for thermal evolution calculations of Swift J174805.3-244637. 


Parameter (unit) / Model 

Standard model 

Shallow heat 

Lower base + shallow heat 

Impurity parameter, Qimp 

1.0 

1.0 

1.0 

Mass-accretion rate during outburst, Mob (g 

1.0 X IQi’^ 

1.0 X lOi’’ 

1.0 X 10^^ 

Outburst duration, tob (yr) 

0.15 

0.15 

0.15 

End of the outburst, to (MJD) 

56166 

56166 

56166 

Atmosphere temperature during outburst, Tob (K) 

2.8 X 10® 

2.8 X 10® 

2.8 X 10® 

Core temperature, Tcore (K) 

1.55 X 10® 

1.55 X 10® 

1.45 X 10® 

Additional shallow heat, Qshallow (MeV) 


1.4 

0.85 

Depth of shallow heat, P/g (g cm~^) 


1.0 X IQi® 

1.0 X IQi® 

xl (dof) 

0.4 (4) 

1.5 (4) 

0.1 (4) 


Note. — The three different models are shown in Fig.|4] See Section lZ3l for details. 


With this lower core temperature, we can set a 2g upper limit of 
Qshaiiow 2.4 MeV nucleon”^, i.e., ~1 MeV nucleon”^ higher 
than obtained for our standard model. Firmly establishing the base 
temperature during the current quiescent phase of Swift J174805.3- 
244637 is thus very important to accurately constrain the amount 
of shallow heating in this neutron star. 


2.4 Empirical decay fits to the temperature curve 

Observationally, the crust-cooling curves of KS 1731-260, XTE 
J1701-462, EXO 0748-676, and IGR J174 80-2446 can be de¬ 
scribed by a (broken) power-law decay (e.g., Cackettet "ZlboiObl: 


scribeg by a (broken) power-law decay (e.g..IC-ackett et alJIzUi^ 
iFridriksson et al.l[201 ll : lOegenaar et alJlioisJ 20141) . whereas 


exponential decay was found to be an adequate, or e ven better, de¬ 
scription for MXB 1659- 29 and MAXI J0556-332 dCackett et all 
l200d : iHoman et al.ll2014h . To compare the crust-cooling curve of 
Swift J 174805.3-244637 with that of these other sources, we there¬ 
fore fitted the temperature evolution to an exponential decay func¬ 
tion of the form y{t) = a + b, and a power-law decay 

of the form y{t) = a (f — fo)~“ + b- Here, a is a normalization 
constant, b a constant offset that represents the quiescent base level 
(i.e., reflecting the core temperature of the neutron star), r the e- 
folding time, a the decay index, and to the start time of the cooling 
curve. 

We performed fits both by fixing the base temperature at the 
value measured before the outburst = 89.7 eV), and by 

leaving it free. The different decay fits are shown in Fig.|5]and the 
resulting fit parameters are given in Table |5] For comparison with 
other sources we also quote the results for a power-law decay fit 
without a constant offset. All fits provide an adequate description 
of the present data; due to the limited number of data points and 
relatively large error bars it is not possible to prefer one statistically 
over the other. An exponential fit favours a base level that is close 
to the temperature measured before the outburst (~ 89.0 eV; red 
dotted curve in Fig. O, whereas a power-law fit tends to a lower 
base level (~ 71.0 eV; black dashed curve in Fig.|5}. 


2.5 A comparison with IGR J17480-2446 

Crust cooling has been detected for six different sources before (Ta- 
ble[TJ. Five of these exhibited a long outburst (> 1 yr), providing the 
conditions for the crust to be significantly heated so that the subse¬ 
quent cooling is observable. However, there is also strong evidence 
for crust cooling in the neutron sta r IGR J17480-2446 following 
an :^2.5 month accret i on outb urst foegenaar & Wiinan dsl l2011bl: 
IPegenaar et al.l[^l Ibl . 120139) . Given the similarly short outburst 
duration, it is therefore interesting to compare Swift J174805.3- 
244637 to this source in particular. 


Table 5. Decay fits to the quiescent temperature curve of Swift J174805.3- 
244637. 


Fit parameter (unit) 

Value 

Exponential decay, base level fixed 

Normalization, a (eV) 

10.1 ±3.3 

Decay time, r (days) 

59.3 ±41.9 

Constant offset, b (eV) 

89.7 fixed 

xl (dof) 

0.39 (3) 

Px 

0.76 

Exponential decay, base level free 

Normalization, a (eV) 

10.6 ±2.9 

Decay time, r (days) 

77.7 ±49.1 

Constant offset, b (eV) 

89.0 ± 1.0 

xl (dof) 

0.20 (2) 

Px 

0.75 

Power-law decay, base level fixed 

Normalization, a (eV) 

75.9 ± 10.2 

Decay index, a 

0.91 ±0.08 

Constant offset, b (eV) 

89.7 fixed 

xl (dof) 

0.70 (3) 

Px 

0.55 

Power-law decay, base level free 

Normalization, a (eV) 

35.8 ±4.9 

Decay index, a 

0.11 ±0.03 

Constant offset, b (eV) 

70.9 ±22.1 

xl (dof) 

0.16(2) 

Px 

0.78 

Power-law decay, no constant offset 

Normalization, a (eV) 

104.1 ±3.7 

Decay index, a 

0.03 ±0.01 

Constant offset, b (eV) 

0 fix 

xl (dof) 

0.20 (3) 

Px 

0.89 


Note. - The start of the cooling curve was set to to = MJD 56166. 


IGR JI7480-2446 is also located in the globular cluster 
Terzan 5 (see Fig. [T]( and therefore has similar Chandra coverage 
as Swift J174805.3-244637. This 11-Hz X- ray pulsar exhibited an 
outbur st in 2010 that lasted for fob — 77 d dPegenaar & Wiinanda 
I 2 OI Ibh . wit h an estimated average accreti on rate of M ~ 2.3 x 
10^^ g s“^ jPegenaar & Wiinandsll20lT^ . Puring the first three 
years in quiescence (2011-2013), the source showed a gradual de¬ 
crease in temperature following a power-law decay with an index 
of —a = 0.47 ± 0.05, and a constant offset of 6 = 77.3 ± 1.0 eV, 
which is slightly higher than the measured pre-outburst tempera¬ 
ture of kT°° = 73.6 ± 1.6 eV dPegenaar et al.ll2013bl, see also 
Fig. |6}. Here, we update the cooling curve of IGR J17480-2446 
using the two new observations of 2013 J uly and 2014 July, whic h 
were obtained after the previous report bv IPegenaar et al.l l l2013bh . 

For the new data we followed the same analysis steps as for 
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Figure 5. Evolution of the neutron star temperature after the 2012 outburst 
of Swift J174805.3-244637 compared to empirical decay fits. The black 
solid (power-law) and red dash-dotted (exponential) fits assume that the 
source settles at its pre-outburst level (2003-2012; grey shaded area). The 
black dashed (power law) and red dotted (exponential) curves instead rep¬ 
resent decay fits with the quiescent base level left as a free parameter. The 
vertical dotted lines indicate the different years in the campaign. 


previous observations of IGR J17480-2446. We fitted the spectral 
data to an NSATMOS model with M = 1.4 M©, J? = 10 km, 
D = 5.5 kpc, and N n = 1.98 x 10^^ cm“^ (for details, see 
IPegenaar et al.l[^13hh . This resulted in neutron star temperatures 
of kT°° = 85.4 ± 1.9 eV (2013 July) and kT°° = 82.7 ± 1.0 eV 
(2014 July). The spectral fits do not require the inclusion of a 
hard emission component, suggesting a contribution of <5% to the 
unabsorbed 0.5-10 keV flux. The updated cooling curve of IGR 
J17480-2446 is shown in Fig. (red data points). The two new 
measurements yielded a temperature similar to that of the preced¬ 
ing observation (2013 February). This is consistent with the slow 
temperature evolution expected at this stage. Indeed, fitting the up¬ 
dated cooling curve to a power-law decay gives —a = 0.45 ±0.04, 
and a constant offset of 6 = 77.3 ± 0.9 eV (shown as the solid red 
curve in Fig. nil. These decay parameters are essenti ally the same as 
obtaine d previously for this source (quoted above; IPegenaar et al.l 
l2013bh . suggesting that in 2013-2014 the source remained on the 
cooling track seen in 2011-2013. 

To provide a direct comparison, we also show the post¬ 
outburst temperature evolution of Swift J174805.3-244637 in 
Fig. |6] (black data points and power-law decay curve). This illus¬ 
trates that shortly after the outburst Swift J174805.3-244637 was 
less hot than IGR J17480-2446, which can plausibly be explained 
by its shorter outburst duration (fob — 2 versus ~2.5 months) 
and lower mass-accretion rate (Mob — lx 10^^ versus ~ 2 x 
10^^ g s“^). Moreover, the temperature curve of Swift J174805.3- 
244637 is flatter, which could be the result of its higher base level 
compared to IGR 117480-2446 ~ 90 versus ~ 74 eV; 

dashed and dotted horizontal lines in Fig. O. The fact that both 
sources show a continuous decay and that the temperature evolution 



Figure 6. A comparison between the two short-duration transients for 
which crust cooling is obseiwed (both located in the globular cluster Terzan 
5). Black triangles are for Swift J174805.3-244637, red filled circles are 
for IGR J17480-2446. Their pre-outburst base levels are indicated by the 
black dashed and red dotted horizontal lines, respectively. The solid curves 
indicate power-law decay fits with a free base level. 

of I GR J1748Q-2446 can b e well described by crust-cooling mod¬ 
els jPegenaar et al .l2013bl) , lends support to the hypothesis that we 
also observe crust cooling in Swift J174805.3-244637. 

Since Swift J174805.3-244637 and IGR 117480-2446 are 
both in the core of Terzan 5 and hence their distances should be 
the same, we can explore how their heat input during outburst and 
the output during the cooling phase compare (with the caveat that 
neutron star parameters also play a role in both the heating and 
the subsequent cooling). To this end we integrated the bolometric 
thermal flux in quiescence from t = 0 till t = 700 d since the 
outburst end (which corresponds to the time of the last available 
data point for Swift J174805.3-244637), after subtracting the con¬ 
stant level due to the core temperature. We then obtain a ratio of 
the fluence in the cooling curve of IGR 117480-2446 compared to 
Swift 1174805.3-244637 of /cool — 5. This is comparable to the 
ratio of the fluences of the outbursts of the two sources (/ob — 3), 
as would be expected within the crustal cooling paradigm. 


3 DISCUSSION 

We reported on Chandra monitoring of the transient neutron star 
LMXB Swift 1174805.3-244637 in the globular cluster Terzan 5, 
following its ~8-week long 2012 discovery outburst. Observations 
were obtained in five different epochs in 2012-2014 and compared 
to seven X-ray spectra obtained prior to the accretion outburst in 
2003-2012. All quiescent spectra are dominated by a soft, thermal 
component, while a non-thermal emission tail contributes ~20-30 
per cent to the 0.5-10 keV unabsorbed flux. Analysis of the pre¬ 
outburst s pectra revealed small variations in this hard power-law 
emission dBahramian et al .l2014h . but we detect no strong variabil¬ 
ity in this spectral component for the post-outburst data. 
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Notably, the fractional contribution of the power-law spec¬ 
tral component to the total 0.5-10 keV flux was consistent with 
being constant in 2012-2014, despite an observed ~35 per cent 
variation in neutron star temperature (see below). This suggests 
that the quiescent thermal and power-law components are phys¬ 
ically linked in Swift J174805.3-244637. Similar results were 
obtained, for instance, for the neutron star LMXBs Cen X-4 


JCackett et al. 2010al: Bemardini et alJ 20131) and MAXI J0556- 

332 jHoman et al. 

2014). whereas in other sources such as 

XTE J170I-462 

EridrikssonetalJl20Ilh and EXO 0748-676 


vealed a clear connection between the two spectral components. 
A tight link between the quiescent thermal and non-thermal emis¬ 
sion could suggest that (at least in some sources) both arise from 
residual accr etion, although this would require fine-tuning accre¬ 
tion models dCackett et al.ll20ld3 : IChakrabartv et al.l[2014f) . Alter¬ 
natively, the power-law tail may be an integral part of the neutron 
star spectrum t hat arises because of the limitations of current atmo¬ 
sphere models jHoman et al. Il20li) . 

Analysis of the seven pre-outburst observations indicated that 
the temperature o f the neutron star in Sw ift J174805.3-244637 re¬ 
mained constant ( iBahramian et al.ll2014h . However, we find that 
~2 weeks after the end of the 2012 outburst the neutron star tem¬ 
perature was significantly elevated compared to that measured be¬ 
fore it became active. The lack of strong variability in the non- 
thermal spectral component implies that there is no obvious indica¬ 
tion that this elevated temperature was due to continued low-level 
accretion. We instead propose that the crust of the neutron star in 
Swift J174805.3-244637 was heated during its 2012 accretion out¬ 
burst and that our Chandra observations sampled the subsequent 
cooling in quiescence. 


3.1 Crust cooling in Swift J174805.3-244637 

If the crust of the neutron star was indeed significantly heated dur¬ 
ing the recent outburst, then our analysis suggests that it cooled 
rapidly, within ~100 d of entering quiescence. Fits to the tem¬ 
perature curve indicate an exponential decay time-scale of ~60- 
80 d, which is considerably shorter than the >400 d cooling time- 
scale of KS 1731-260 and MXB 1659-29 (both these sources 
continue to cool many ye ars after their outburst ended; e.g., 
ICackett et alj |200A l2010bl) . However, considering the typically 
large errors on the characteristic decay times (see e.g., Tabled, 
the cooling time of Swift J174805.3-244637 is not significantly 
different from the ~ 150-200 d decays observed for EXO 0748- 
676, IGR J17480-2446, and M AXI J0556-332 jPegenaar et al.l 
l2013bl , [ 2 OI 4 I ; iHoman et alj|2014t) . On the other hand, the data of 
Swift J174805.3-244637 may also be consistent with continued 
cooling following a power-law decay, which would suggest a post¬ 
outburst base level ~20 per cent lower than observed prior to the 
2012 accretion phase. This could potentially be explained as a dif¬ 
ferent amount o f FI/He left on the surface of the neutron star after 
its last outburst ( iBrown et ^l2002l : iMedin & Cummindl2015L see 
also Section[T2](. 

A direct comparison with crust cooling observed after a 10- 
week accretion outburst of IGR J17480-2446 suggests that the 
more rapid decay and flatter temperature curve of Swift J 174805.3- 
244637 may naturally be explained by its higher base tempera¬ 
ture. As a result, the temperature profile in the crust at the end of 
the outburst is much less steep than for a longer outburst and/or 
a lower core te mperature, so that the thermal relaxation is muc h 
quicker (see e.g.. lBrown & Cumming||200'^ ; |Page & Reddvll2012h . 


This may also be the reason that there is no clear evidence for (long¬ 
term) crust cooling in Aql X-1 l lBrown et alJll99d ; iRutledge et al.l 
l2002ah . Its outburst properties are quite similar to Swift J 174805.3- 
244637 (on average fob — 70 d. Lx ~ 3 x 10®® ergs“®; 
ICampana et al.l2013ll . but its q uiescent temperature is considerably 
higher, kT°° ~110 eV (e.g.. ICackett et all 1201 ih . Neutron stars 
with low pre-outburst quiescent temperatures may therefore be the 
best targets for future studies, even if their outburst is short. 

It is of note that rapid temperature decays wer e seen at the end 
of the (short) outbursts of e.g., XTE J1709-267 dPe genaar et H] 
12013^, IGR J17494-3030 dArmas Padilla et alJl2013ll FI and Aql 
X-1 Jcampana et^l2014l) . The characteristic time-scale in these 
sources is of the order of a few days, which is comparable to the 
temperature chan ges seen during quiescen t accretion flares in e.g., 
XTE 11701^62 dFridriksson et al.ll20lih . and MAXI J0556-332 
dHoman et al .I 2 OI 4 I) . The similar, short decay times in these sources 
could perhaps be caused by draining (part of) the accretion disc on 
a viscous time-scale, rather than cooling of the neutron star crust 
dFridriksson^tay|20]T[ jArmas Padilla^talj 201_3|;JCampana_et_al.l 


|2014 Homan et al. 2014 ). However, Medin & Cummin j ~i 2015h 

suggested that rapid cooling within ~1 week after the end of an 
accretion outburst can be explained as cooling of the outermost 
(fluid) layers of the neutron star called the ‘ocean’. These authors 
show that this model ca n explain the rapid tem perature evolution 
seen in XTE J1709-267 dPegenaar et al.ll20r3^ . However, the de¬ 
cay in the very-faint LMXB IGR J17494-3030 was strikingly sim¬ 
ilar to that of XTE J1709-267, despite a factor of ~100 difference 
in outburst accretion rate. A more prominent difference in tempera¬ 
ture evolution might therefore be e xpected, perhaps arguing ag ainst 
the cooling envelope interpretation dArmas Padilla et al.l2013 ). The 
observed cooling time-scale for Swift J174805.3-244637, on the 
other hand, is much longer than the rapid temperature decays seen 
in these sources, and is fully consistent with crust cooling. 


3.2 Shallow crustal heating 


The occurrence of shallow heating in the crusts of neutron stars, not 
accounted for by standard nuclear reactions, remains a puzzle. An 
additional heat source of ~l-2 MeV is required to properly model 
the crust-cooling curves of MXB 1659-29, IGR J17480-2446, and 
EXO 0748-676 (ferown & Cummin3 l2009l: IPegenaar et akll^l Ibl ; 

m 


IPage & Reddvll2013l : Pegenaar et all2014 ) , whereas a lower value 

of ~ 0.1-0.5 MeV is suffi c ient fo r KS 1731-260 and XT E J1701- 
462 dBro wn & Cu mmind d2009ll : IPage & Reddvl d2013h . but see 
iTurlione et alj I ^15l) ) □ Furthermore, a number of neutron stars 
that display superbursts and frequency drifts in mHz quasi-periodic 
oscillations also appea r to require additional sha l low heating (e.g., 
Cumminget^ 200^ fee^t all I 2 OO 8 L l2009l : lAltamirano et al.l 
2008Ll2012l:lLinaresetalj2012l). 


It is currently unknown whether shallow heating occurs for all 
neutron stars and, if so, whether the magnitude of the extra energy 
release should always be the same. The small number of sources 
and limited constraints on shallow heating in the current sample 
of crust-cooling curves do not allow us to determine whether the 


^ IGR 117494-3030 is a very faint LMXB, not to be confused with the 
crust-cooling source IGR J17480-2446 in Terzan 5. 

^ We note that the crust-cooling curve of MAXI J0556-332 has not been 
modelled yet, but the very high temperatures observed for this source 
iHoman et alj2014h will likely also require the addition of substantial shal¬ 
low heating. 
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lack/requirement of extra heat is related to source-specific prop¬ 
erties such as e.g., mass, spin, magnetic field strength or the type 
of companion star. One mechanism that could perhaps account for 
additional heat fl ux in shallow crustal layers is compositionally- 
driven convection jMedin & Cummin j201 iLboiShlPegenaar et"^ 
l2013ah . This process strongly depends on the composition of the 
liquid ocean (the remains of nuclear burning on the surface of the 
neutron star). Therefore, it may give rise to differences in shal¬ 
low heating for different sources, and a single source may also 
experi ence different levels of h eating from one outburst to an¬ 
other jMedin & Cummin3l2015h . This could potentially be tested 
with crust-cooling observations, e.g., by investigating if the tem¬ 
perature curves of different sources can all be explained with the 
same amount of shallow heating, or by comparing the crust-cooling 
curves of a single source after different outbursts. 

The crust-cooling curve of Swift J174805.3-244637 can be 
adequately modelled using standard physics input, without the need 
to invoke an additional source of shallow heat. However, a heat 
source up to ~1.4 MeV nucleon”^ is still be compatible with the 
present data (at 2a confidence), which is not unlike that found for 
MXB 1659-29, IGRJ17480-2446, and EXO 0748-676. 

An important caveat is that our conclusions on shallow heat¬ 
ing in Swift J174805.3-244637 are based on the assumption that 
the neutron star has (nearly) settled at a base temperature similar 
to that observed prior to the 2012 outburst (i.e., implying that the 
crust has fully cooled). However, the temperature of the neutron 
star does not necessarily have to level off at exactly the same value 
as observed before the outburst. The amount of unburned H/He on 
the surface of a neutron star determines the heat flux that flows from 
the crust. Theoretically, for the same core temperature, the observed 
surface temperature could therefore differ by a factor of a few afte r 
different outbursts jBrown et al.ll2002l : iMedin & Cummin3l2015h . 
although there is no strong evidence that this effect is indeed ob¬ 
servable. Empirical (power law) decay fits allow for the possibil¬ 
ity of continued crust cooling in Swift J174805.3-244637, which 
would change the derived constraints on the shallow heating. Eur- 
ther observations of Terzan 5 would allow us to confirm/reject that 
the crust continues to cool, and hence to solidify the constraints on 
shallow heating in this neutron star. 
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